function fu = f(u)

fu = u*u/2;

end